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Abstract 

This paper describes the effort to provide a preliminary capability analysis and a neural- 
network based adaptive flight control system for the JPL-led BEES aircraft project. 

The BEES flyer was envisioned to be a small, autonomous platform with sensing 
and control systems mimicking those of biological systems for the purpose of scientific 
exploration on the surface of Mars. The platform is physically tightly constrained by the 
necessity of efficient packing within rockets for the trip to Mars. Given the physical 
constraints, the system is not an ideal configuration for aerodynamics or stability and 
control. The objectives of this effort are to evaluate the aerodynamics characteristics of the 
existing design, to make recommendations as to potential improvements and to provide a 
control system that stabilizes the existing aircraft for nominal flight and damaged conditions. 
Towards this several questions are raised and analyses are presented to arrive at answers to 
some of the questions raised. 

CART3D, a high-fidelity inviscid analysis package for conceptual and preliminary 
aerodynamic design, was used to compute a parametric set of solutions over the expected 
flight domain. Stability and control derivatives were extracted from the database and 
integrated with the neural flight control system. The Integrated Vehicle Modeling 
Environment (IVME) was also utilized for estimating aircraft geometric, inertial, and 
aerodynamic characteristics. 

A generic neural flight control system is used to provide adaptive control without the 
requirement for extensive gain scheduling or explicit system identification. The neural flight 
control system uses reference models to specify desired handling qualities in the roll, pitch, 
and yaw axes, and incorporates both pre-trained and on-line learning neural networks in the 
inverse model portion of the controller. Results are presented for the BEES aircraft in the 
subsonic regime for terrestrial and Martian environments. 
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1. Introduction 


Fast-moving airplanes designed to navigate Mars the way insects control their courses on Earth might one 
day travel to places where exploration rovers and landers cannot reach. Mars is covered with intriguing 
canyons, cliffs, and craters, most of them too remote or hazardous for conventional exploration vehicles. 
Exploration Air Vehicles (EAVs) could carry cameras and other sensors to these spots. Building a Mars EAV 
will provide a new way for NASA to conduct exploration. In addition to exploration of hard to reach spaces, 
EAVs could relay information between rovers already exploring the planet's surface or monitor a spacecraft 
while it touches down. 

The JPL-led Mars EAV project is named BEES, for Bio-inspired Engineering of Exploration Systems, and 
has been in the works for a few years [1,2]. The plane is called a biomorphic flyer because it captures some 
of nature's successful flight designs — in this case, unique strategies of navigation, hazard avoidance and 
flight— and combines them with conventional designs. The BEES technology is based on the fact that insects 
like the honeybee and the dragonfly possess a brain with less than 0.01 percent the neurons found in a human 
brain, yet with simple strategies they navigate in the world just fine. 

Putting a plane on Mars poses several challenges. For one, the atmosphere on Mars is so thin that flying an 
aircraft near the surface is difficult, like trying to control a plane 100,000 feet above the Earth's surface. And 
because Mars is mostly unexplored, an unpiloted aircraft would need to have the ability to sense and respond 
to what it encounters. 

The work reported in this paper examines software based on an adaptive neural network to control the BEES 
EAV’s attitude in the midst of uncertainties and failures. 

2 . Test Article 

The current BEES flyer is shown in Figure 1. The aircraft is a delta-type wing flyer with twin vertical tails 
for stabilization. Two elevons are utilized for control of the aircraft. Information for the 5kg biomorphic 
flyer was utilized in this effort. The values are referenced below: 


Table 1: Test Article (5 kg, 02 Biomorphic Flyer) Details 


Item 




Dry Mass 

2670g 

Landing Speed 

=55km/h 

Fuel 

340g 

Range 

55000m 

Payload 

1990g 

Endurance 

>20min 

Takeoff weight 

5000g 

Wing Span 

1.240m 

Engine. ... 

698g 

.ivi/iv-f 


Propeller 

58g 

Leading edge sweep 
back 

50.8° 

Engine Type 

OS61FX 

Aspect Ratio 

2.206 

Propeller Size 

28.0cmx20cm 

Taper Ratio 

0.19 

Top Speed 

=165km/h 

Wing Area 

0.701m' 
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Figure 1 : The BEES Platform 


3. .Modeling Using CART3D 


CART3D is a high-fidelity inviscid analysis package for conceptual and preliminary' aerodynamic design [3. 
4], It allows users to perform automated CFD analysis on complex geometry. The package includes utilities 
for geometry import, surface modeling and intersection, mesh generation and flow simulation. Currently the 
package includes two flow simulation codes: Tiger, and flowCart. The package is highly automated so that 
geometry acquisition and mesh generation can usually be performed within a matter of minutes on most 
modem UNIX workstations or PC's. 


Geometry comes into the package in the form of surface triangulations which can be generated from 
many popular CAD packages, from legacy surface triangulations or from structured surface grids. C.ART3D 
uses adaptively refined Cartesian grids to discretize the space surrounding a geometry and cuts the geometry 
out of the set of "cut-cells" which actually intersect the boundary. The flow solvers run on single and multi- 
processor computing platforms. 

A Pro-E solid model of the BEES aircraft was provided by Australian National University (ANU) as 
the basis for this effort. Substantial efforts were involved in removing features (such as the engine and holes 
for the actuation devices) and in making the geometry “water-tight. The engine cowling was modified to 
provide a complete seal against the body. The slot for the propeller was not faired over for this effort. Figure 
2 shows a sample of the CAD geometry' as generated by the Pro-E solid modeler. A view of a representative 
grid of the final geometry is provided in Figure 3 and demonstrates the locally adaptive Cartesian nature of 
the grids utilized by CART3D. Looking carefully at Figure 3, it can be seen that the geometry contains 
significant surface waviness up to about 20 % of chord which could affect some features of the aerodynamic 
field. 
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Figure 2. BEES CAD Model - First Cut 



Figure 3. Cartesian Grid Section on Computational Geometry for Fully Deflected Elevons. 


3.1. Results Using CART3D 

The analysis using CART3D generated 4704 cases for the following combinations of flight variables 
Mach # = 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.75, 0.8 
Angle-of-attaek = -10, -5, 0. 5, 10, 15 degrees 
Side Slip Angle = 0,-10 degrees 

Elevon Deflection = 49 different cases (each elevon at -10, -5, 0, 5, 10, 15 and 20 degrees); 
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Surface pressure data was integrated to provide forces and moments on the delta wing, stabilizers and control 
surfaces of the BEES fiver. These were used to analyze and answer the fundamental question: Will the BEES 
Flyer fly at subsonic speeds in Mars? 

Figures 4 through 7 document several of the aerodynamic curves of interest. 

Given the plots, one can compute the velocity needed to fly level, non-accelerating flight at different positive 
angles-of-attack. 


Table 2: Mach # and Thrust Required for Level Flight 


Given: 

Mass of Flyer (m) — n kg 

Gravity on Mars (g M ) = 1/3 of Earth gravity (g E ) 

p = Density on Mars at "sea level" = 0.015 kg/m ’ 

a» = Speed of sound on Mars =268 m/s 

S = surface area = 0.7 m~ 

Cl = Coefficient of lift 
Cd = Coefficient of drag 
L = Lifting force 
D = Drag 

W = Weight of craft on Mars = mg M 
T = Engine thrust 
a = Angle of attack 
V = Forward airspeed 


Assumptions: 

No wind or turbulence 

Thrust required is to maintain level flight only 


Angle-of-attack 

(deg) 

Velocity (m/s) 

Mach # 

Thrust Required (N) 

5 

122.18 

0.46 

0.71 

10 

83.88 

0.31 

1.55 

15 

69.67 

0.26 

2.57 1 
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CD Vs Alpha For Varying Mach # 



alpha (deg) 


Figure 4. C D Versus Alpha for Varying Mach # 


CL Vs Alpha For Varying Mach # 



alpha (deg) 

Figure 5. C L Versus Alpha for Varying Mach # 
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Cm Vs Alpha For Varying Mach # 



-10 -5 0 5 10 

alpha (deg) 

Figure 6. C m Versus Alpha for Varying Mach # 


CY Vs Alpha For Varying Mach #@ beta=-10 deg 


cn cn -*r in co r-. 




By examining the data presented, one can conclude the following: 

a. The BEES flyer can fly at low angles of attack in the subsonic regime with no needed modifications. 

b. The drag rise for Mach > 0.7 is quite severe even at moderate angles of attack. This problem could 
be alleviated by either avoiding the transonic regime, or preferably, by modifying the design. A 
small change in the airfoil could result in a significant reduction in drag in the transonic regime. 

c. The engine casing and the hole behind it distorts the flow significantly. This needsrto be redesigned. 
The initial estimate of control authority seems sufficient for nominal flight conditions. 

d. Sensor locations have to be chosen to minimize the effects of the flow. 

3.2. Extensions to the model to Include Rotational Effects 

The CART3D solutions generated for this study are steady in time and thus don't directly permit the 
computation of stability and control derivatives with respect to the pitch, roll and yaw rates of the aircraft. 
Unsteady solutions can be generated but take significantly more computational resources. An approximation 
to these derivatives can be formed by assuming that the effect of rotation is dominated by the change in flow 
angle relative to the body as one moves away from the center of gravity. For the pitch and roll axes, the 
effect would then be dominated by the derivative with respect to angle of attack: 
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And for the yaw axis, it is assumed that the primary effect will be the sideslip of the vertical surfaces relative 
to the freestream: 
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4. Intelligent Flight Control System 

Figure 8 presents the implementation of a Level 1 Intelligent Controller for the BEES test bed developed in 
this study. It should be noted that Level 0 is non-adaptive whereas Level 1 is adaptive. The generic neural 
flight control system is based on the augmented model inversion architecture developed by Rysdyk and 
Calise [11]. This direct adaptive tracking controller integrates feedback linearization theory with both pre- 
tiained and on-line learning neural networks. Details of the aichitecture are given below. 

Reference Models: The autopilot commands roll and pitch rates and aerodynamic lateral accelerations 
through elevon and rudder inputs. These commands are then transformed into body-axis rate commands, 
which also include turn coordination, level turn compensation, and yaw-dampening terms. First-order 
reference models are used to filter these commands in order to shape desired handing qualities. 

P + I Error Controller (Level 0): Errors in roll rate, pitch rate, and yaw rate responses can be caused by 
inaccuracies in aerodynamic estimates and model inversion. Unidentified damage or failures can also 
introduce additional errors. In order to achieve a rate-command-attitude-hold (RCAH) system, a 
proportional-integral (PI) error controller is used to correct for errors detected from roll rate, pitch rate, and 
yaw rate (p, q, r) feedback. 

Learning Neural Network (Level 1): The on-line learning neural networks work in conjunction with the error 
controller. By recognizing patterns in the behavior of the error, the neural networks can learn to remove 
biases through control augmentation co mm ands. These commands prevent the integrators from having to 
windup to remove error biases. By allowing integrators to operate at nominal levels, the neural networks 
enable the controller to provide consistent handling qualities. The learning neural networks not only help 
control the nominal system, but also provide an additional potential for adapting to changes in aircraft 
dynamics due to control surface failures or airframe damage. 

Dynamic Inversion/Aero Generation: The dynamic inversion element converts the summed response 
commands into virtual control surface commands. Dynamic inversion is based upon feedback linearization 
theory. No gain-scheduling is required, since gains are functions of aerodynamic stability and control 
derivative estimates and sensor feedback. Several methods are available to accomplish approximate model 
definition: simple linear model methods, nonlinear tables or using pre -trained neural networks (non-adapting) 
to provide estimates of aerodynamic stability and control characteristics. The model is then inverted to solve 
for the necessary control surface commands. In our work, a Levenberg-Marquardt (LM) multi-layer 
perceptron is used to provide dynamic estimates for model inversion. The LM network is pre-trained with 
stability and control derivative data generated by a Rapid Aircraft Modeler, and vortex-lattice code [5]. 
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Figure 8. Intelligent Flight Control Architecture 

5. Results and Discussion 

Results are documented for the following cases: 

• Mars and Earth scenarios with open-loop control and IFC system 

• Control redundancy options with open-loop control and IFC system 

o Two elevons and no rudder 

O JTVJU1 CiCVOilU Willi I1VJ iUUUCi 

• Failure cases with open-loop control and IFC system 

o Failed elevons (all combinations) 
o Longitudinal and lateral CG shift 

Figure 9 presents roll and pitch doublet data for the BEES flyer in both Earth and Mars atmospheres at 0.5 
Mach number for both an open loop control and the neural flight control system. Yaw data is not presented 
because without rudder panels, the yaw response of the aircraft cannot be directly controlled. These doublets 
were generated using the perturbation equation option of the FLTz code. As can be seen in the figure, the 
aircraft response in Earth’s atmosphere is generally well behaved for both open-loop and IFC control 
systems. In the Mars atmosphere, undamped oscillations greatly affect the command tracking for both pitch 
and roll response in the open loop case. For the IFC controller, the aircraft response is reasonably consistent 
between the two environments. Figure 10 presents similar data for a 4 elevon configuration with no failures. 
The response characteristics are similar to those of the two elevon case. 

The first of the failure cases presented in Figure 1 1 is that of a two elevon configuration with the left elevon 
jammed. With one of the two available control surfaces failed, independent control of pitch and roll rates is 
no longer possible. Apparent stability on Earth (left plots) depends strongly on starting in a trim condition. 
The immediate divergence in roll when a pitch input is applied in the Mars open loop case may not be 
immediately obvious. IFC fights to retain control, but the best it can achieve is a trade between pitch and roll 
rates. In Figure 12, a four elevon case with the left outboard elevon failed is considered. Splitting the control 
surfaces gives the error controller the opportunity to regain and maintain control after a single actuator 
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failure. Note the early termination of the simulation in the open loop cases as the craft tumbles through 
attitudes outside the simulation domain. 

In Figure 13, a longitudinal CG shift is imposed on the two eleven configuration. The longitudinal CG shift 
changes the pitch response of the system; in this case, the CG moved toward the rear of the craft by 5% of 
the MAC. The open loop response on Earth shows increased pitch response, which is damped back down to 
reference levels by the error controller. The open loop response on Mars shows increased instability. Note 
that despite the large initial response to the doublet, IFC brings this condition back under control. In Figure 
14, the same longitudinal CG shift is applied to the four elevon configuration. Splitting the control surfaces 
neither improves nor degrades system response to longitudinal CG shift. 

The two elevon configuration with a lateral CG shift failure is shown in Figure 15. The lateral CG 
shift induces a huge coupling from the longitudinal control into the lateral response axis, as can be seen by 
the -initial response on Earth of the open loop system to a pitch input. The error controller brings this 
response back under control. The Martian conditions exacerbate the instability of the system. The error 
controller works in both environments to counteract the coupling, bringing the rate back under control and 
largely returning the attitude (not shown) to expected values. Early termination of the open loop simulations 
was due to the tumbling of the craft through attitudes outside the domain of the simulation; the transient at 41 
seconds in the IFC earth case was a controlled descent through the terrain. The four elevon lateral CG shift 
case is shown in Figure 16. 

Splitting the control surfaces neither improves nor degrades system response to a lateral CG shift. 
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Figure 11. Aircraft Response Doublets for Two Elevon Configuration Left Elevon Failed. 
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Figure 13. Aircraft Response Doublets for Two Elevon Configuration. Longitudinal CG Shift 
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6. Conclusions 


Our analysis using a low-fidelity math modeling tool in conjunction with the CART3D high fidelity CFD 
modeler showed that even in an undamaged configuration, stability and control of this particular 
BEES configuration will be difficult in the Martian atmosphere without the aid of an error controller. In 
addition, any uncertainties will certainly require an adaptive controller. Conclusions are presented for control 
and design issues. 

• Control issues 

o Need of EFC for Mars application was established from both stability and controllability 
points of view. 

o Additional control effectors (analytical control redundancy), such as splitting the elevators 
and adding rudders could improve control effectiveness and resilience in the face of damage. 

• Design Issues 

o CART3D provided initial verification of flyer's abilities. Use of a different wing planform to 
reduce drag rise at higher Mach numbers and angles-of-attack should be investigated. 

7. Recommendations 

Our recommendations for future studies related to design and control are presented next. First we present 
highlights of our recommendations and then present details of design trade-off issues and the outer-loop 
guidance and control issues. 

• Design Issues 

o The drag rise for Mach > 0.7 is quite severe even at moderate angles-of-attack. This problem 
could be alleviated either by avoiding the transonic regime, or preferably by modifying the 
design. 

o Sensor locations have to be chosen to minimize the effects of the flow, 
o Flyer design and the controller design should be approached simultaneously towards the goal 
of achieving the mission needs. 

• Control issues 

o Outer-loop guidance and control issues are important and were not examined in this study, 
o Sensor noise/failure issues need to be examined. 

o Hardware implementation of the IFC design using Field Programmable Gate Arrays (FPGA) 
need to be investigated. 
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